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Abstract 

The low- momentum effective interaction Vi ow k has been formulated in the three-dimensional 
momentum-helicity representation as a function of the magnitude of momentum vectors and the 
angle between them. As an application, AV18 potential has been used in the model space of 
Lee-Suzuki method and it has been shown that the low-momentum effective interaction, Vi ow k 
reproduces the same two-body observables obtained by the bare potential Vjvjv- 
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I. INTRODUCTION 



Several methods have been developed to derive the energy independent low-momentum 
effective interaction such as the renormalization group (RG) and the model space tech- 
niques. These approaches are mainly based on a partial wave (PW) decomposition and 
the details have been given in references Recently Bogner et al. have developed a 

low-momentum effective interaction which describes the two-nucleon system at low energy 
successfully. This effective interaction is independent of the potential models as the cutoff 



is lowered to A = 2.1 fm~ l [lo|, 111] 



During the past years, the three-dimensional (3D) approach has been developed for few- 
body bound and scattering problems (^ol . In this approach momentum helicity basis 
states have been used for the representation of the nuclear forces. The motivation for devel- 
oping this approach is to introduce a direct solution of the integral equations avoiding the 
very involved angular momentum algebra occurring for the permutations, transformations 
and especially for the three-body forces. Conceptually the 3D formalism considers all partial 
wave channels automatically. 

Recently we have developed a 3D formalism for construction of low-momentum effective 
interaction neglecting the spin and isospin degrees of freedom [27 1. Considering the spin 
and isospin is a major additional task, which we intend to present in this article. Our aim 
is to extend the low-momentum effective interaction directly in a spin-isospin dependent 
3D approach and to formulate the low-momentum effective interaction in the momentum- 
helicity representation. 

This article is structured in the following way: In section[II]the model space of Lee-Suzuki 
method has been used to derive the energy independent model space effective interaction in 
the 3D momentum-helicity representation. In section HTT1 the reduced forms of the equations 
have been displayed by choosing suitable coordinate systems. Section [TV] describes the 
numerical calculations of the low-momentum effective interaction Vi ow k in the model space 
Lee-Suzuki method by using the AV18 potential. Finally a summary and an outlook have 
been provided in section [Vj 
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II. LEE-SUZUKI METHOD IN THE 3D MOMENTUM-HELICITY REPRESEN- 
TATION 



The Lee-Suzuki method has been applied to the free space nucleon-nucleon problem in 
the 3D momentum-helicity representation and the low-momentum effective interaction Vi ow k 
has been obtained as a function of the momentum vectors. In the model space methods the 
projection operators onto the physically important low-energy model space, the P space, 
and the high-energy complement, the Q space, have been introduced on momentum-helicity 
basis state |k; kS'A; t)™ as: 

P = I d\t\k]\tS\]ty a j ™(k;kSX;t\, |k| < A, 

SXwt J 

Q = [ dk\k;kSX;ty a ^ ™{k;kSX;t\, \k\ > A, (1) 

where A is a momentum cutoff which divides the Hilbert space into the low and high momen- 
tum states. The antisymmetrized momentum-helicity basis state which is parity eigenstate 



is given by 



25): 



\k;kS\;t)™ = -±=(l-P 12 )\k;kSX) n \t) 



= ^(l-vA-n^SX^lt), (2) 

Here S is the total spin, A is the spin projection along relative momentum of two nucleons, 
t is the total isospin and \t) = \trrit) is the total isospin state of the two nucleons. m t 
is the isospin projection along its quantization axis which reveals the total electric charge 
of system. For simplicity m t is suppressed since electric charge is conserved. In Eq. (J2]) 
P12 is the permutation operator which exchanges the two nucleons labels in all spaces i.e. 
momentum, spin and isospin spaces, and |k; kS'A),,- is parity eigenstate which is given by: 

|k;LSA>^ = -J=(l + Vn P n )\k;kSX), (3) 

where P n is parity operator, rj n = ±1 are the parity eigenvalues and |k; kS'A) is momentum- 
helicity state. The normalization of the momentum-helicity basis state is given by 251 ] : 

^(k^k's'A'^'l^kSA; t)* a = (l-vA-) S+t )St'tS v „, v Js>s 

x {<5(k' - k)8 yx + ^(-)^(k' + k)5y- x }, 

(4) 



and the completeness relation of this state is defined by: 

[ dk l k ; kSA; t)™ ^ ™(k; kSX; t\ = 1. (5) 

It is clear that the projection operators P and Q satisfy the following relations: 

P + Q = 1, 

PQ = QP = 0, 

P 2 = P, 

Q 2 = Q, (6) 

and they act on the full-space two-body problem states as: 

WW™ = ^i^r"- 

Q\^sx; t r = <o\**fr>J m , (7) 

where I^VkSA-t) 7 ™ an< ^ l^k-icSAt) 7 ™ denote the states of the full and model spaces respectively 
and lo is an operator which transforms the states of the P-space to the Q-space. The non- 
hermitian low-momentum effective potential in the model space that reproduces the model 
space component of the wave function from the full-space wave function is given by Qj-js]: 

V lowk = PV NN (P + QuP), (8) 

where Vjvjv denotes the bare two-body interaction. By using the integral form of the pro- 
jection operators P and Q, the low-momentum effective interaction Vi ow k can be written in 
the 3D representation on momentum-helicity basis state as: 

™( V ';f>'SX;t\V lowk \ V -i>S\-ty a = m (p';-p'SX';t\V NN \p-pSX;ty a 

+ r dqq2 I dci na (p';t>'SX';t\V NN \q;qSX";ty a 

x na (q;qSX";t\oj\p;pSX;ty a , (9) 

where we use the property that Vnn conserves parity, spin and isospin. By defining matrix 
elements of each arbitrary operator A as: 

^ 5t;A ' A (p', p) = ™(p'; p'SX'; t\A\p; pSX; t)™. (10) 



4 



Eq. (Q can be rewritten as: 

T/ 7r5t,A'A/ / \ T/ 7r5t,A'A/ / \ . 1 f , 2 f i~ r rirSt,\'\" / i \ irSt,\"\/ \ /i i \ 

V low k {P,P) = V NN ' (P,P) + 7^ / d( lQ d <\ V NN (P,q)w ' (q,p),(ll) 

where p and q are the momentum vectors in the P and Q spaces respectively. As we know, 
the total spin states of two nucleons are singlet (S = 0) and triplet (S = 1) states. For the 
singlet case, Eq. (TTTj) . there is only one equation which reads as: 

V lowk (P,P) = Vnn (p,P) + ^J^ dQQ J dqV NN ' (p , q) U) ' (q,p), (12) 

however, for the triplet case there are three equations for each initial helicity A = — 1,0, 1 
as: 

T /7rlt.A'A/ / \ T r7rlt,A'A/ / \ . 1 f j 2 f t« r T /7rlt,A'l/ / \ 7rlt,lA/ \ 

Kw k (P'P) = V nn (P-PJ + ^y^ dqq j dq{V NN > (p , q) u ' (q, p) 

. T r7rlt,A'0/ / \ 7rlt,0A/„ \ i T/7rlt,A'— 1/ / \ 7rlt,— 1A/ \i ;,n\ 

+ *W (p.q)^ (q.p) + tW (p,q)w ' (q,p)}. (13) 

Symmetry relations for the potential matrix elements in the momentum-helicity representa- 
tion are denoted by |25J: 



t rTtSt,— 

NN 


A'A 


y,p) 


= ^(-) 5 ^' VA 


(-p', 


p 


T/ 7rSt,A 
NN 




y,p) 


= ^(-« A ' A 


&>',- 


p 


NN 




y,p) 


= ^ A ' A (- P ',- 


p). 





(14) 

A corresponding symmetry relations are also valid for cu-matrix elements. By applying this 
symmetry relations for matrix elements of V^n and u, Eq. ffT3]) can be reduced to: 

V£? X (P',P) = V^ X '\ V \ V ) + \f™ dqq 2 J dq{2^ A,1 (p',q)a;^ 1A (q,p) 

+ ^ A ' (p',q)a;^ 0A (q,p)}. (15) 

Hence for the triplet case one needs only two equations for A' = 1, for each A. To calculate 
the low-momentum effective potential Vi ow k we need to determine cu- matrix elements. The 
key aspect of the Lee-Suzuki method is the determination of the u operator defined by the 



following equation 



28|: 



Q\A&x*) m = Q"P\A&>#)' m ' ( 16 ) 
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Applying """(q; qSA'; t\ to the left hand side and using the integral form of projection oper- 
ators P and Q, this equation can be rewritten as: 

™(q;qSA';t|^ k; k5A;i)™ = f I ft™ t\u>\p-,pS\ u ;t) m 

x ™(p;pS\' , ;t\* k . Mx ., t y a . (17) 

We introduce the completeness relation for the scattering wave function in the model space 
on momentum-helicity basis state: 

\ n J J 

= (1 - ^(-)^){5( p ' - p)5 yx + Vn (-) s 5(p' + p)^_ A }, (18) 



where we use the following notations for the matrix elements of the wave functions: 

^,yA (p) s -( P ;p5A';^ k;fcsA;t )-, 

^,AV (p) s ^.^Jp^A';*}™ (19) 
After implementing Eq. ([TBI . Eq. (ITTj) can be written as: 

^st,vA (q)P) = y A d ^ y d k^,A'A" (q) ^,A"A (p); (20) 

where \l/£ St ' A A ( P ) and \P k s *' A A (q) are the wave function components of the P and Q spaces 
of the full-space, respectively. These can be written in the form of the half-on-shell (HOS) 
two-body T-matrix as: 

K St ' X "(P) = (1 - vA~) S+t mP - k)5v A + vA-) S S(P + k)5v-A + - , 2 J }, 

2 k z — p z + is 

(21) 

^A'A (q) = T^q^g) (22) 



The HOS two-body T-matrix can be obtained from the Lippmann-Schwinger (LS 



equation 



in the 3D momentum-helicity representation, which is given by (h = = 1) 25]: 

T^ A ' A (k',k,^) = y;r A (k',k ) + ix: / ** NN k^ + ie ' 

(23) 
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The equation for the singlet case is: 



-i7r0t,00/i,' 1, 7„2\ _ T/ 7r0t,00/,/ ,^ , 1 / j\JI NN l K ' K I 1 l K ' K > K ) 



7«(k', k, = v^(k', dk" -^ fc ;_fc,, 2+ ,r ' (24) 

and for the triplet case we have two coupled equations for A' = 1, for each A as: 

1 f ^ ll V^ x '\k'X)T^(k l ',k,k 2 ) 



T^ VA (k',M 2 ) = ^> A ' A (k',k) + I J dk" 



k 2 - k" 2 + is 



+ 4 dk k*-k>» + ie • (25) 



Eqs. f fT8l) and ( 1201 for the singlet case can be written as: 

r-A 
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±j\kk 2 Jdk ^i 0t ' 00 (p')% 0tfi0 (p) = (i - vA-Y)W - p) + vAp' + p)}, (26) 



^•°°(q, P) = \£ dk k 2 J dk vfjf ' 00 (q) ^ M ' 00 (p). (27) 

It is clear that the same symmetry relations as we mentioned for potential matrix elements 
are valid for matrix elements of \l/ and Thus by applying the symmetry relations for 
triplet case, we write two equations for A' = 1, for each A respectively: 

i j* ^ dkk 2 J dk{2% lt ' x '\p')% lt ' lx (p) + ^ u - A '°(p')^ u ' 0A ( p ) } 

= (1 + V7T (-f){6(p' - p)6 yx - r]J(p' + p)<5 A '-a}, (28) 

^"' A ' A (q,p) = \ f dkk2 j rfk{2^ A,1 (q)^ 1A (p)+ ^"-^(^^"^(p)}. 

(29) 



III. CHOOSING SUITABLE COORDINATE SYSTEMS FOR NUMERICAL CAL- 
CULATIONS 

In this section suitable coordinate systems are chosen in order to write the Eqs. fill I) . (JTSJ) 
and ( l20i) in reduced forms for numerical calculations. The azimuthal behavior of potential 
and two-body T-matrix elements for special case where the vector k is along z axis arc given 

by m 

T^ A ' A (k',fcz) = e 4V T" 5t ' A ' A (fc', £;,£'), 

Vf N ^'\k',kz) = e^'vf N ^'\k',k,x'), (30) 
7 



where x' — k ■ k. Inserting these relations into the LS equation one can obtain [25l |: 

rpnSt,X'X(ul t. T l t,2\ _ ynSt,X'X,y , /\ 

+1 E jf ik " k " £ ^" ,V >'*''^^J^'' A ' fc, '' : '"''' :2 > . (31) 

where x" = k • k and (//'-integration is carried out independently by defining: 

uf^'^'ik', k", x', x") = dip" e -*W-V') ^' A ' A "(k', k"). (32) 

Jo 

The integrand is periodical with respect to p" . Thus we can set ip' = and Eq. (I3TI) for the 
singlet and the triplet cases can be written as: 

T 7r(tt,00 (jj u I h 2\ _ v nOt,00 (h , , A 

J. yrv ^ tv j Jj j tv J — NN \ 5 5 / 



1 



dk"k"j^dx"^ — v ' 7; , ' ' ' 7 , issi 



T 



■Ti-lt.A'A n„' 7„ ™/ J.2 

Jb ; 

OO /*1 

2 



+1 / rf^'A;" / dr" 
o </-i 
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x")T n0tfi0 {k", 


k. 


x", /c 2 ) 




k 2 


- k" 2 + ie 




j 












vl U > X '\k>,k" 




,x")T* u > lx (k" 


, h 


x", fc 2 ) 




k 2 


- k" 2 + ie 






ul u ' X '°(k', k" 


,x> 


,x")T* u > 0X (k" 




' ) 



-1 /»oo /» 1 

+ W ^ ^^"j,;..- I (34) 



J-l 

irSt,X'XnJ 1, J,2\ 



As we know the T-matrix elements T ' (k , k, ) are not the solution of the LS equa- 
tion. The solution of the LS equation would be T nSt,xx (k',k,x',k 2 ), which are T-matrix 
elements in momentum-helicity basis with initial momentum in the z-direction and without 
its azimuthal dependence. Therefore T'^'W, k, k 2 ) can be connected to the solution of 

n 

the LS equation as follows [261 ] : 

T^ x ' x (k\k,k 2 ) = ^ N =~ s f NX ' [X)dNx[X) T^ x ' x (k',k,y,k 2 ), (35) 

d \'x(y) 

where: 

y = k ■ k = x'x + \J\ — x' 2 \J\ — x 2 cos(<p' — <p), (36) 



and df, y (x) are rotation matrices 



'A'AV 

o;-matrix elements: 



29(. Now we consider the azimuthal behavior of ^ and 



co^ x ' x (qi,p) = e- iX 'W St > X ' x (q,p,x), 

% st ' x ' x (qz) = e- M ^l st ' X ' x (q,x), (37) 



S 



where the vector q is along the z axis and x = p • q. The Eq. ( 120|) can be reduced to 



S* ,xx (q t p,x) = t dkk2 f (Lr " n S! - yy '«h-'-") k^"\p,x" ,). 

4 A „ Jo J-i 



where x" = k • q and: 

V k (q,x) 



«St,X>X», T^^"(q,k,x",k 2 ) 



k 2 — q 2 

Integration over tp" can be performed independently by defining: 

¥ k s } x " x (p,x",x) = / d/ e ^(^)ff^(p), 
Eq. (1381) for the singlet and the triplet cases can be written as: 

W *<*«W, *) = \[ dk £ dx" x") x", x), 



u* u > x ' x (q,p,x) = I J dkk 2 £dx"{2*l U > X '\q,x") ¥ k ^ ix (p, x" , x) 
+ ^ k u ' X '°(q,x") ¥ k u y °\p,x",x)}. 

Multiplying both sides of Eq. (|T8l) by e~ tX " & and integrating over <p' yields: 

W t dkk2 t dx"^; x ' x "( P ',x',x") ^(p/.j) 

4 A „ Jo J-l 

= (1 - vA~) S+t ) 5{p' - p){S(x' - x)5 yx + V7T (-) S+X '"5(x' + x)8 x ^ x }, 



where: 



/•2tt 

Jo 



and: 



n st,xx (po = (i-^(-) 5+t )w-fc)[5(x'-x")w-^')^A" 

EL-5 ^'"^ 4a»(^) r^V, k, y, k 2 ) 

2d s yx „(y) k 2 -p' 2 + ie 



Eq. (jUJ) for the singlet and the triplet cases can be written as: 

I 



^ /" ji (2 / 7 // ,T.7r(M,CI0/ / / //\ .T r vr0t,00/ // \ 

-y dA;£; J dx V kfi ' (p , a; , x ) * fc)0 ' {p,x ,x) 



(1 - r^R'W - p){6{x' -x)+ rjj(x' + x)}, 



A f 1 

dkk 2 d x" {2 ^fc^/'// 1 (p', x', x") ^l^y} x (p, x" , x) 
o J-i 



x{S(x' - x)5 yx - ^{-) x '"5{x' + x)5y. x }. (47) 

Finally by considering the vector p along z axis and using azimuthal behavior of potential 
and cu-matrix elements, Eq. ( Till can be rewritten as: 



Vio S w f\p'^) = V^ x '\p>,p,x') 

I p oo pi 

+ lE/ ^%f t ' A ' A ''(p',g,x',x'')^ A '' A (g,p,x''), (48) 

where x' = p' • p and x" = q • p. This equation for the singlet and the triplet cases can be 
written as: 

1 poo pi 

+ i/ A M'^tf t q,afy)u,^{q, P y), (49) 

V lowk (P,P, X ) = V NN {P,P,X) 
i poo pi 

+ 4 J d 1<l 2 J i dx " i 2 ^'V, *") ^"' 1A (g,P, x") 

+ ^ u ' X '°(p',q,x',x")uj^ \ q ,p,x")}. (50) 

IV. DISCUSSION AND NUMERICAL RESULTS 

We have chosen AV18 phenomenological potential for our calculations. This potential is 
fitted to pp as well as np data below 350 MeV laboratory energy. In addition the AV18 
potential is fitted also to low-energy nn scattering parameters and deuteron properties. With 
this interaction in the first step we have calculated two-body T-matrix by solving the LS 
Eqs. f l33|) and (jMj) for the singlet and the triplet cases respectively [25(. In the next step we 
have calculated ^j^m A (p, x",x) as an inverse of ^f^m X (p',x',x") for each A'" from Eqs. 
(|46p and (|47p by using the LU decomposition method. In the numerical calculations we have 



used the Lapack library 30[], for to solve a system of linear equations for the calculation of 
^fcV" 1 X (p,x",x) and the two-body T-matrix elements. Then by solving Eqs. (jUl) and (142]) 
we have obtained u 7rSt ' x ' x (q,p, x) and finally we have inserted the a;-matrix elements into 
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Eqs. (HHD and (ISTij) to obtain the low-momentum effective interaction V^^ x A (p',p, x') for 

the singlet and the triplet cases respectively. 

In numerical calculations we have used the Gaussian quadrature grid points to discrete 

the momentum and the angle variables. The integration interval for the P and Q spaces are 

covered by two different hyperbolic and linear mappings of the Gauss-Legendre points from 

the interval [-1,+1] to the intervals [0, a] and [a, b]\ 

l+x u_b- a b+a 

k - 1 71 2yZ' k-——x-\ — . (51) 

A ^ A a> X Z Z 

The typical values for a and b are 10 fmT 1 and 150 fwT 1 , respectively. As we mentioned 
in the introduction section we have used the value of 2.1 fm~ l for the cutoff A in our 
calculations. The ^/'-integration within an interval [0, 2n] has been rewritten within an 
interval [0, ~] as shown in the following notation: 

/ = dip"f(cos((f' - ip")) e imi ^'-^ = dip" f (cos p") e imv " 
Jo Jo 

= r dp" {/ '(cos tp") e imip " + /(- cos p") e im ^" +7T) } 
Jo 

dip" { f (cos ip")(e im *" + e im{2 ^" ] ) + /(- cos(^")(e fm( ^" ) + e im{ ^" ] }. (52) 



o 

The second equality has been justified by the periodicity of the integrand within 2n. Thus 
the number of integration points for polar angle has been reduced. 

In our calculations we have chosen sixty grid points for the momentum variables in the 
interval [0, a], and twenty two grid points for the momentum variables in the interval [a, b}. 
Also forty and ten grid points for the spherical and the polar angle variables have been used 
respectively. The solutions of the integral Eqs. (!33|) and AMI) require a one-dimensional 



3 II ] for its accuracy and 



interpolation. We have used the cubic hermitian splines of Ref. 
high computational speed. 

In Fig. [T]we have shown the total two-body cross section for np scattering by using Vi ow k 
and Vnn- Also in Fig. [2] we have compared the calculated 1 S'o and 3 Si-wave phase shifts 
from Vi ow k and Vnn- For calculation of the phase shifts we have used the relation between 



the PW and 3D representations of the on-shell two-body T-matrix 



26|: 



T S ^' l (p, P y) = \ V^^ 1 C(l'Sj; OA') C(lSj; OA) 



2? , 

J X'X 



X 



J dx' di x ,(x') T* st > yx (p,p, x',p 2 ). (53) 
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1000 



E 




FIG. 1: The total cross section from the low-momentum effective potential Viowk (circles) and the 
bare potential Vnn (solid line) as a function of kinetic energy in the lab frame. 




FIG. 2: The 1 S and 3 5i -wave phase shifts from the low-momentum effective potential V/ oli > & 
(circles) and the bare potential Vnn (solid line) as a function of kinetic energy in the lab frame. 

The results are in good agreement with high accuracy. In Figs. [3}{6] the calculated low- 
momentum effective potential V^\ x x (p, p, x') and the bare potential V^' x x (p, p, x') as well 
as differences between them have been shown as a function of the momentum variable p and 
the angle variable x. 

As a test of our calculations we have compared the obtained results for the low-momentum 
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P [fm _1 ] 



-1 

(a) 



P [fm" 1 ] 



-1 

(b) 



FIG. 3: (a) The comparison of the low-momentum effective potential Vi ow k (solid lines) with the 
bare potential Vnn (dashed lines) and (b) differences between them, for S = and t = 0. 




FIG. 4: (a) The comparison of the low-momentum effective potential Vi ow k (solid lines) with the 
bare potential Vnn (dashed lines) and (b) differences between them, for S = and t = 1. 

effective potential and the uo operator in the 3D and the PW approaches. As a first step 
we have calculated the uj operator and the low-momentum effective potential in the PW 
approach directly. We have then obtained the PW projection of the u operator and the low- 
momentum effective potential from their corresponding 3D representation by the following 
relations: 
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FIG. 5: (a) The comparison of the low-momentum effective potential Vi ow k (solid lines) with the 
bare potential Vnn (dashed lines) and (b) differences between them, for 5 = 1 and t = 1. 
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FIG. 6: (a) The comparison of the low-momentum effective potential Vi ow k (solid lines) with the 
bare potential Vnn (dashed lines) and (b) differences between them, for S = 1 and t = 0. 
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FIG. 7: The comparison of ui{q,p){p 2 — q 2 ) for the 1 So-wave calculated in the PW approach (dotes) 
and in the 3D approach (solid lines). 



vSX'toP) = \ v ^2j + E CVSj ' 0A '> c " Sj; 0A> 



A' A 



x J^dx' d{ x ,(x>) V^Y\p,p,x'), (54) 
* W W) = \ ^ T ^+ J C(VSj; OX') C{lSr, OA) 



J A'A 
1 

7r5t,A'A 



x J dx d{ x ,(x) u wbt ' x \q,p,x). (55) 

The obtained results for the channel 1 ^ in the 3D and PW approaches have been given 
in Fig. [7] and Fig. [HJ The agreement between the two approaches is quite satisfactory. 
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FIG. 8: The comparison of the low-momentum effective potential V lowk for the 1 S'o-wave calculated 
in the PW approach (dotes) and in the 3D approach (solid lines). 

V. SUMMARY AND OUTLOOK 

In this article the 3D formulation of the model space Lee-Suzuki method based on 
momentum-helicity representation has been presented on moment um-helicity basis states. 
The low-momentum effective interaction Vi ow k has been derived as a function of the magni- 
tude of momenta and the angle between them for the singlet and the triplet cases respectively 
without using the partial wave decomposition. The calculated two-body observables from 
the low-momentum effective interaction and the bare interaction have been presented. In 
addition, a comparison between the calculated Vi ow k from the PW and the 3D approach has 
been demonstrated as a test of our calculations. 

The advantage of our formulation in the 3D representation in comparison with the 
PW representation is that we have calculated the low-momentum effective interaction by 
considering all partial waves automatically. The implementation of the obtained 3D low- 
momentum effective interaction in the few-body bound and scattering calculations is a major 
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task that can be considered for the future investigations. 
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